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We consider the influence of the dark energy dynamics at the onset of cosmic acceleration on 
the Cosmic Microwave Background (CMB) bispectrum, through the weak lensing effect induced 
by structure formation. We study the line of sight behavior of the contribution to the bispectrum 
signal at a given angular multipole I: we show that it is non-zero in a narrow interval centered at a 
redshift z satisfying the relation l/r(z) ~ kNhiz), where the wavenumber corresponds to the scale 
entering the non-linear phase, and r is the cosmological comoving distance. The relevant redshift 
interval is in the range 0.1 z 2 for multipoles 1000 ^ I ^ 100; the signal amplitude, reflecting 

the perturbation dynamics, is a function of the cosmological expansion rate at those epochs, probing 
the dark energy equation of state redshift dependence independently on its present value. 

We provide a worked example by considering tracking inverse power law and SUGRA Quintessence 
scenarios, having sensibly different redshift dynamics and respecting all the present observational 
Oh' constraints. For scenarios having the same present equation of state, we find that the effect described 

above induces a projection feature which makes the bispectra shifted by several tens of multipoles, 
about 10 times more than the corresponding effect on the ordinary CMB angular power spectrum. 
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INTRODUCTION 



The recent cosmological observations indicate that the Universe is nearly geometrically flat, filled with structures 
which grew out of a primordial linear spectrum of nearly Gaussian and scale invariant perturbations; the scalar contri- 
y—i ,; bution, i.e. from the perturbations in the energy density, appears to dominate over the tensor one from cosmological 
00 ' gravitational waves. About 5% of the critical energy density is made of baryons, while the remaining dark part is 
supposed to interact at most weakly with the baryons themselves, since we observe it only indirectly, i.e. through 
' its gravitational effects. The dark component appears to be 30% pressureless, like in Cold Dark Matter (CDM) 
I scenarios, dominating the gravitational potentials perturbations which host visible structures like galaxies or clusters. 
• The remaining 70% should be in some sort of vacuum energy, with negative pressure acting as a repulsive gravity, 
Oh' and responsible for a late time cosmic acceleration era. 

Q . Despite of the remarkable convergence of the most important cosmological parameters, the picture is clearly far 
$H ' from being satisfactory: without a better insight into the nature of the dark cosmological component, we cannot 
c/3 claim to have a satisfactory physical understanding of cosmology. In particular, the evidence for a non-zero vacuum 
energy driving cosmic acceleration today (hereafter dark energy) is a clear sign that the study of the dark cosmological 
^ ] component could have unexpected and most important consequences on cosmology and more generally on fundamental 
physics. The evidence for the dark energy is rather robust, since it comes from two independent datasets. The first one 
is represented by the observations of distant standard candles, namely type la supernovae (hereafter SNIa), directly 
. probing the cosmological expansion rate at low redshifts l] ; these observations indicate that a 70% fraction of the 
■ ■ ' critical energy density today is in the dark energy. The same result is supported, independently, by the observations of 
the Cosmic Microwave Background (CMB) anisotropics and Large Scale Structure (LSS). The abundance of clustered 
matter is known to be at the 30% level of the critical density by observing the nearby Universe jj]; on the other 
hand, the scale subtended in the sky by the acoustic oscillations in the CMB anisotropy angular power spectrum 
clearly support a flat cosmological geometry (see the recent results by the Wilkinson Microwave Anisotropy Probe, 
map.gsfc.nasa.gov, hereafter WMAP, Q and references therein). Thus, CMB and LSS together indicate again a 
70% of critical energy density in the vacuum. The only way to break this evidence would be invoking a low value of 
the Hubble constant, which is excluded by the Hubble Space Telescope (HST) data 0]. 

In cosmology, a constant vacuum energy density is provided by the Cosmological Constant. While a dark matter 
component has well established support from the theories beyond the standard model of particle physics, a Cosmo- 
logical Constant at the observed level raises two fundamental questions. The first is the fine-tuning: the energy scale 
is about 120 orders of magnitude less of the Planck energy density which is supposed to be the unification scale of all 
forces in the early Universe. The second question is the coincidence: not only the vacuum energy scale is small with 
respect to any fundamental scale, but it is such that it is comparable to the matter energy density today. 

Already several years ago 0, 01 it was realized that the evidence for a Cosmological Constant would have raised 
these crucial questions. To address them, the Cosmological Constant was minimally extended to a dynamical and 
inhomogeneous scalar field, evolving slowly enough to yield an almost constant vacuum energy out of its potential. 
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providing cosmic acceleration. The latter evidence caused a renewed interest in these models f^,^- In particular, it 
was demonstrated how the dynamics of this component, under suitable potential energy, can originate attractors in 
the trajectory space, named tracking solutions, capable to reach the present dark energy density starting for a in the 
remote very early Universe, thus alleviating, at least classically, the problem of fine-tuning mentioned above 
The scalar field playing the role of the dark energy was named Quintessence. This description allowed to investigate 
the relation of the dark energy with the other cosmological components: while an explicit coupling with baryons is 
severely constrained by observations JJJ , the relation of the Quintessence with the gravitational ^3 ^^'^ dark matter 
[13| sectors of the fundamental Lagrangian, as well as the phenomenology arising from generalized kinetic energy 
terms Q have been extensivel y st udied. The Quintessence field has been proposed to be related to super-symmetry 
and super-gravity theories (see |l5l | and , respectively, and references therein) but we still lack a reliable candidate 
from fundamental physics. 

Generally speaking, despite of the large number and variety of ideas proposed to explain the evidence for cosmic 
acceleration, going from quantum spacetime microphysics to exotic particle physics (see [T3 and references therein), no 
valid solution to the coincidence problem has been provided yet. It is likely that a breakthrough in the understanding 
of the dark energy could come from the future high precision data on the redshift dependence cosmological acceleration. 
Indeed, almost all the candidates for the dark energy predict a characteristic behavior of the equation of state, w, 
as a function of the redshift z (see and references therein). A good example of this occurrence, which we will 
use extensively in this work, are the different Quintessence potentials proposed so far, causing relevant differences in 
w{z) at redshift relevant for cosmic acceleration, even if its present value is close to the Cosmological Constant case. 
Wo = -1 _19]. 

The recent CMB observations, combined with the LSS, HST and SNIa data, allow to constrain the effective dark 
energy equation of state, assumed to be constant in time, to be Weff< —0.78 at 2<t confidence level imposing w > —I 
prior and Weff = —0.98 ± 0.12 dropping the previous prior (see [2C| and references therein). That will be certainly 
improved by future observations, but the real challenge in order to discriminate between different models is to gain 
insight into the dark energy equation of state redshift dependence. Together with the observations of distant SNIa, 
it is likely that the study of cosmological structure formation will be crucial for investigating the onset of cosmic 
acceleration, since the redshift intervals relevant for the two processes overlap. The best way to look at the structure 
formation process is just to study the weak gravitational lensing deflection induced by forming structures along the line 
of sight of the background light (see (2lj for reviews) . The shear induced by weak lensing on background galaxies has 
been observed by several groups, and the results agree impressively between different tele scop es and data reduction 
techniques [23; already, several dark energy weak lensing observables have been proposed [2J. New powerful probes 
to map the weak lensing shear over large sky areas are being designed to operate in the next decade with ground-based 
and space surveys (see 24] and references therein). On the CMB side, the weak lensing effect on the total intensity 
and polarization CMB anisotropics can be predicted on the general basis of cosmological perturbation theory . The 
effect is generally non-Gaussian, and relevant at the arcminutes angular scale; it is likely that the low noise, arcminute 
resolution imaging expected by the Planck satellite (see astro.estec.esa.nl/SA-general/Projects/Plajick/) and 
future CMB missions (see CMBpol spacescience.nasa.gov/missions/concepts.htm) will have a crucial impor- 
tance in weak lensing studies. 

The weak lensing effect on the CMB is the subject of the present work. We focus on the third order statistics, 
where the background CMB light, if well described by a Gaussian process at last scattering as it appears from the 
recent CMB observations j26rj . yields a zero contribution within cosmic variance. As we already mentioned, a suitable 
harmonic approach to the weak lensing effects on the CMB total intensity and polarization has been carried out [2^ , 
out of a general description of the CMB anisotropics [23|. Further works focused specifically on the effects on the 
CMB bispectrum 28], and in particular on the benefits of having bispectrum data in addition to the CMB power 
spectrum, allowing more precise estimates of the main cosmological parameters, even breaking the degeneracies when 
the dark energy is included into the analysis |29| . In this work we investigate the sensitivity of the bispectrum to the 
redshift dependence of the dark energy equation of state w{z). 

In Section^ we review the basic features of the Quintessence dynamics in tracking trajectories, which we will use 
later to illustrate the phenomenology we study. In Section Hill we review the CMB third order statistics, described 
through the bispectrum, and how it is affected through the weak gravitational lensing induced by structure formation. 
In Section llVl we study the bispectrum line of sight integral describing how it is sensitive to the cosmological expansion 
rate redshift dependence, and in Section we provide a first quantification of the effect by considering the leading 
Quintessence scenarios. Finally, in Section IVII we make the concluding remarks. 
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II. QUINTESSENCE COSMOLOGY 



The dynamics of the Quintessence scalar field cj) in cosmology follows the Klein-Gordon equation, describing back- 
ground and linear perturbation evolution: 

D(p + V' = , S{a(j>) + V"S(j) = , (f ) 

where the □ is the Dalambertian operator in a Friedmann Robertson Walker (FRW) cosmology, and 6 represents 
linear fluctuation, while V{(p) is the Quintessence potential and V' — dV/dd. 

We consider here inverse power law (hereafter RP, first considered in Q) and SUGRA ^f^ potentials, defined 
respectively as 

Vrp = , VsuGBA = ■ exp (4^G02) . (2) 

The exponential in the SUGRA model induces relevant changes in the field dynamics with respect to the RP case, as 
we explain now. As we mentioned in the Introduction, the tracking trajectories in the Quintessence scenarios alleviate 
the fine-tuning in the early Universe, since by means of them one is no longer forced to assume that the dark energy 
is vanishingly small with respect to radiation in the early Universe [SIEIalSE|: the present field value (/)o, of the 
order of the Planck mass Mpianck, is reached from a wide set of initial conditions and 0^, with the only important 
condition that (pi ^ Mpianck- Until the latter regime is satisfied both in the RP as well as the SUGRA cases, the 
tracking trajectories are well defined until the Quintessence is subdominating, yielding an almost constant equation 
of state obeying the simple relation 

PDE 2 

^-TT—^ (3) 

PDE Z + a 

where pde = {dcj) / dt)"^ / 2 — V and pde = Pde + 2V are the field pressure and energy density, respectively. In the 
RP scenario, if we require the present equation of state to be in agreement with the current constraints y, say 
— 1 ^ Weff ^ — 0.8, the exponent is in the range ^ a ^ 0.5, yielding a shallow potential shape. At the present, 

the tracking regime is abandoned because the dark energy is no longer a subdominant component, but the shallow 
potential shape makes the present equation of state not far from the tracking one in Eq. differing typically at the 
10% level. The SUGRA exponential correction flattens sensibly the potential shape at ^ ~ Afpianck, i-e. at the end of 
the tracking trajectory. This means that a given equation of state at present is obtained for values of a sensibly larger 
than in the RP case, meaning that the dark energy dynamics and thus the cosmological expansion rate as a function 
of the redshift are generally much different in the two cases. This aspect is phenomenologically illustrated in Fig. ^ 
featuring the redshift evolution of the equation of state in SUGRA and RP tracking trajectories yielding at present 
Wo = —0.9, from the present to the tracking regime; we consider a cosmological scenario which is close to the current 
best fit 0, i.e. a flat Friedmann Robertson Walker (FRW) cosmology with Hubble parameter h = 0.7, baryon, CDM 
and dark energy fraction density flth^ = 0.022, Hcdm + ~ 0.3, floE = 0.7, respectively, three species of massless 
neutrinos and a scale invariant scalar perturbation spectrum with no gravitational waves. The behavior of w(z) in 
the figure is markely different in the two Quintessence scenarios considered, refiecting the two values of a required to 
have wq; specifically, a = 0.34 and 1.76 in the RP and SUGRA cases, respectively. We will also consider an effective 
dark energy model with constant w — —0.9 (hereafter DEe//) obtained simply by giving that equation of state to the 
Cosmological Constant term in the FRW background cosmological evolution. That is physically incorrect because it 
assumes no spatial fluctuations in the dark energy; the latter are present indeed, and play a relevant role even in the 
minimally coupled Quintessence scenarios on super-horizon scales, as it can be seen by looking at the large scale CMB 
anisotropies [33; in addition, the dark energy spatial fluctuations are present on all scales if the field interacts with 
other cosmological components 0,0,^3- The curves in Fig. as well as all the perturbation spectra quantities 
in the next Section, have been numerically computed by making use of a dark energy oriented version (3Lj of the 
CMBFAST code 

The most relevant quantities for the following discussion are the distance and the dark energy density behavior as 
a function of the redshift z. The first is given by 

/ ^ c r dz' 
r{z)^— / , , (4) 



^"•^^ ^JnMo{l + zf + pDE{z')/pco 



where we neglected the radiation contribution and pco is the critical cosmological density at present; pde{z) is the 
redshift dependent background dark energy density, given by 

,\ + w{z') 



Pde {z) = pD_Eoexp 



dz' 
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FIG. 1: Equation of state redshift dependence for the Quintessence models described in the text, for the RP model (dotted 
line) and the SUGRA (dashed line). Note that both have too at the present. The solid line is a DEe// model with a constant 
equation of state. 

as it can be verified by simply solving the energy conservation equation. In the RP and SUGRA scenarios described 
above, t{z) and posiz) have different behavior as a consequence of the different behavior of w{z) shown in Fig. ^ 

We conclude this Section by recalling briefly the main dark energy effects of the CMB angular power spectrum ,3^ . 
As w gets larger than —1, the distance in Eq. gets reduced. This causes a shift of the CMB angular power spectrum 
toward smaller multipoles, i.e. large angular scales, affecting in particular the acoustic peaks location. Moreover, for 
w > —1 the dark energy domination era begins earlier that in for the Cosmological Constant, since increases 
with redshifts. This boosts the gravitational potential dynamics at low redshift, injecting dynamics on large angular 
scales, where however the Cosmic Variance greatly limits the impact of the effect. Thus the projection represents the 
main dark energy effect on the CMB angular power spectrum. However, it should be noted that a limiting factor for 
the CMB power spectrum as a probe of the dark energy is the fact that the CMB path performs a sort of redshift 
average of the dark energy properties along their line of sight, making the CMB angular power spectrum mainly 
sensitive to the redshift average of the dark energy equation of state in the interval in which it is relevant, say z ^ 2, 

quite low with respect to the origin of the CMB at z ~ 1000; the latter statement could be violated if the dark energy 
is not subdominant at last scattering "SSf]. Nevertheless ^Jt should be noted that a great improvement in that respect 
is expected by the next generation CMB satellites (see 01 for the Planck case). 

On the other hand, the injection of power in the CMB anisotropics through the forming structures is potentially 
a stronger dark energy probe, since the relevant redshift interval is close to the onset of cosmic acceleration, as we 
pointed out already. This process is likely to alter the statistical properties of the CMB last scattered radiation . 
In particular, if the latter is Gaussian as it appears from the recent CMB observations |2^, it is convenient to design 
CMB observables looking where the noise due to the last scattered signal is lower, for example in the third order 
statistics. In the following we focus on this aspect. 
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III. WEAK LENSING AND CMB BISPECTRUM 



In this Section we write down the equations describing the weak Icnsing effect induced by structure formation on 
the CMB bispectrum. We decompose the CMB anisotropy in a direction fi into two terms: 



e (n) ~ Qiss {n + d) + Qisw (n) 



(6) 



the first term represents the CMB anisotropy at decouphng [Iss means last scattering surface), which was last 
scattered on a direction n + a and gravitationally lensed to our line of sight fi; the second term is the CMB anisotrOT 
contribution from the formingstructures along the line of sight, the Integrated Sachs- Wolfe effect (hereafter ISW [22 
also known as Rees-Sciama jS^])- We do not consider here the contributions coming from the Sunyaev-Zel'dovich 
effect [3^, which yields different power at different frequencies, and the Ostriker-Vishniac effect [33 • To the first 
order in a, Eq. © becomes 



e 



(n) ~ Qiss (n) + Qisw (") + VOi^s (n) 



a 



and the ISW term can be expressed as 



e 



ISW 



in) 



/d 
dr—-9{r],hr) , 



(7) 



(8) 



where 77 is the cosmological conformal time, r the comoving distance along the line of sight. 



[(1 



2*^ 



and 



4" is the cosmo- 
m dx^] , if the 



logical gravitational potential, defined as the fluctuation in the metric ds^ 
anisotropic stress is zero [2^ [23 ■ 

Note that in general Eq. (jS)) describes the contribution from linear and non-linear density fluctuations, with the 
only assumption that the gravitational potential fluctuations they induce are linear. The effect of gravitational lensing 
on the CMB is contained into the last term of Eq. {T)), indeed, a is the gradient of the lensing potential defined as 
the projection of gravitational potential along the line of sight (see plllSq '): 



dr 



ris 



-4' (r, fir) 



(9) 



Note that the ISW and lensing effects are correlated since they both arise from ^I*; this correlation induces a non- 
Gaussian feature on the CMB pattern and leads to a non vanishing bispectrum signal 25, ^8, ,29,, 39J . 

The gravitational potential is built out of the density fluctuations through the Poisson equation. To describe 
the non-linear part of the power spectrum, we adopt a semi-analytic prescription [i^ which evaluate the non-linear 
power out of a remapping and rescaling of the linear one, as we describe below; note that, as we show explicitly in 
the next Section, our results do not depend on the particular recipe we adopted for treating the non-linear density 
perturbations, which have been described in several ways (see for example |4l| and references therein), but mainly on 
the perturbation dynamics in linear regime. The non-linear (NL) and linear (L) density fluctuation power spectra 
are related as 



PNLik,z) = G 



ATTk^PL{k,z) 
P'^a^iz)^ ■ 



PL{kL,z) , 



(10) 



where the linear one is evaluated at the wavelength 



kL 



[l + 47rfc3PwL(fc,z)]'/' 
the linear power spectrum is defined as usual as 

PL{k,z) = Ak"T^{k,z) , 



(11) 



(12) 



where A is the primordial perturbation amplitude, n is the spectral index, and T{k, z) is the matter trans- 
fer function. Note that Eq. H1U|) with Hll|) represents an implicit equation which has to be solved 
for Pt^L(k,z) for any wavenumber k at any epoch z. In the relation H1U|) G is a function describ- 
ing the rise of the power spectrum at the onset of non-linearity; specifically, a suitable fit is G{x, z) = 
[1 + ln(l + a;/2)] [l + 0.02a;4 + 1.08 • IQ-'^x^ / g{z)] / (1 + 2.1 • IQ-^x^^/^), with (3 = 0.83 and / = .go|wo|i-3l"'°l-°-^s 
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FIG. 2; D factor as function of redshift for our three models: solid line is DEe//, dotted line is RP and dashed line is SUGRA. 
The three curves, almost superposed, do not show significant differences between the three models. 



|40j. We are interested in the case in which the dark energy equation of state changes with time. The latter modifies 
the background dynamics and therefore the perturbation growth rate, defined as usual as 

g(z)=go(l + .) ^(^^Q^Q) ■ (13) 

Thus 'w{z) enters in the Unear perturbation growth rate as well as in the normalization trg as a function of the 
redshift. As we mentioned in the previous Section, all the quantities above, i.e. A, T{k,z) and g{z), as well as 
(Ts{z) are computed numerically in the Quintessence cosmologies of interest here [sij : the perturbation amplitude A 
is determined normalizing the corresponding CMB power spectrum to the WMAP data. 

The CMB bispectrum is built out of the third order statistics in the harmonic domain, and it is defined as 

where aim are the coefficients of the expansion of the total intensity CMB fluctuations into spherical harmonics; 
we choose to work with dimensionless quantities when treating the CMB: in particular the aim are dimensionless 
corresponding to the harmonic expansion of AT/Tq where Tq is the present CMB thermodynamical temperature. 
If the statistics of the signal is rotationally invariant, it is convenient to work with the so called angle-averaged 
bispectrum psf: 

Bi,i,i, = Y ( '0 5,7?^'"^ , (15) 

iii2<3 ^ \ mi 7712 ma / ^ 

mim2m3 ^ ^ 

where the parenthesis are the Wigner 3 J symbols. Expanding in spherical harmonics Eq. Q, the expansion coefficients 
[2^|3^ can be written as 

ISW I / ,\T?x+m'+m" 



aim — aim ~^ "-Im ^ / , \ ) ^-m+rn'+m" 



I'l" 
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FIG. 3: Integrand for SUGRA model with too = —0.9 as function of redshift for different multipoles; solid line is linear regime, 
dotted line is non linear regime. The graphics are qualitative the same for DEe// and RP, the integrand is dimensionless. 



where aj™'* are the coefficients of the harmonic expansion of Eq. 0; in the equation above G^mlm^imz the Gaunt 
integral defined in appendix^ After some algebra (see the appendix in |29|| for a detailed derivation), the bispectrum 
can be expressed by exploiting the expression above as 



(2h + 1) (2^2 + 1) (2/3 + 1) / h h h 
Att 



(17) 



where 5P. indicates the permutation over li,l2 and I3, Cf'^" is the power spectrum of primordial CMB anisotropics, 
and 



Q(o^((aLr)*«Lr) 



2 / dz \'''\ \^ 
Jq r [ziss) r-i [z) 



dPm (k, z) 



dz 



(18) 



where P$ (fc, z) is the gravitational potential power spectrum, related to the density power spectrum by the relation 

P-q,{k,z) — (I^^Mo)^ ("^) P {k, z) {1 + z)'^ , and ziss is the redshift of last scattering surface; Q (l) is the most 
relevant quantity here, describing how the forming structures along the line of sight induce the lensing on CMB 
photons, expressed as the statistical expectation of the correlation between the ISW and lensing effects |2^l39ll4^. 
In particular, the expression above has been used to evaluate the bispectrum dependence on the most important 
cosmological parameters, including an effectively constant dark energy equation of state, and the benefits of the 
bispectrum data on the estimation of the cosmological parameters themselves j2flj . In the next Section we study how 
the contribution to Q(l) is distributed along the line of sight, and how it depends on the dark energy equation of state 
redshift dependence. 
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FIG. 4: Ratio between non linear and linear mass power spectrum for DEe// (solid line), RP (dotted line) and SUGRA (dashed 
line) with wo = —0.9 as function of wavenumber k for two redshift. 



IV. THE LINE OF SIGHT DISTRIBUTION OF THE BISPECTRUM SIGNAL 

The integrand in Eq. p8|l is built out of a geometrical factor D made of distances and a derivative factor F based 
on the cosmological fluctuations, this two factors are defined as 



r [zissjr* [z) 



dP^ (k, z) 
dz 



(19) 



In this Section we describe and study their relevance for probing the dark energy dynamics. 

The factor D is not significantly sensitive to the dark energy redshift dynamics. The reason is that the dependence 
of on z gets redshift averaged through the distance integral 10}; the latter gets the dominant contribution from the 
lowest redshifts, say z ^ 0.1, where all cosmologies converge to the same behavior. Thus, the behavior of w at higher 

redshifts yields a negligible contribution to the distances, even if the differences are relevant as in Fig. ^ In Fig. [21 
we plot D for the three cosmologies defined in Section^ showing explicitly the poor sensitivity of D to the w(z) 
behavior. 

The most important contribution in the integral H18|l comes from the gravitational potential fluctuation derivative 
term F; we now study its behavior along the line of sight. For a given multipole Z, and coming from high redshift 
z, dPqr/dz is probed from large to small scales A(z) = 27rr(z)/Z. Let us study the asymptotic regimes z — > oo and 
z — > 0. In the first case, for multipoles of the order 10^, at very high z the scale A is outside the horizon; in this 
regime, if the Universe is matter or radiation dominated the potential is constant, so that F vanishes. In the second 
case, X —>■ 0, F vanishes again simply because the power spectrum derivative is probed at infinite wavenumbers. 
Thus the relevant bispectrum signal must come from some region in between these two asymptotic regimes. This is 
evident in Fig. |31 where we plot the integrand of Q{1) for different multipoles in the SUGRA scenario described in 
Section m as a function of the redshift. As z comes down from infinity, the scale A decreases so that eventually it 
meets the growing comoving scale corresponding to the horizon, 1/aH; from that point A(z) becomes a sub-horizon 
comoving cosmological scale. In this regime, the gravitational potential decreases first, because of the free-streaming 
of the underlying matter density fluctuations. This feature is visible in Fig. 13 as the tail of the curve at dP-^/dz > 
coming from high z. As z gets smaller, eventually A(z) matches the scale which is entering the non-linear phase at 
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the epoch z\ at this point the derivative changes sign, dPq,/dz < 0, since the power is now increasing with time, i.e. 
as z decreases. Thus, the negative peaks in Fig. O map the non-hnear region of the spectrum. As z is reduced to 
zero, the power vanishes approaching infinite wavenumbers as described above. Thus, the F term yields a non-zero 
contribution to Q{1) in Eq. (|18|l only in a redshift interval around the epoch corresponding roughly to the rcdshift 
Zbs satisfying the relation 



where the left hand side is purely geometric, while kNL represents the scale entering the non-linear regime at Zbs- 

Now, the effect we just described turns out to be most interesting in at least three aspects. First, as it can be seen 
in Fig. 13 for multipoles between 1000 and 100 the solution to Eq. H20|l is in a redshift interval extremely interesting 
for the onset of cosmic acceleration, say between 0.1 and 2, respectively. Second, the peak cuts out the present epoch: 
this means that the bispectrum integral probes the cosmological perturbation dynamics at the relevant redshift for 
cosmic acceleration, independently on the present. Third, for a fixed multipole the amplitude of the signal reflects 
the perturbation dynamics at the epoch Zbs, which in turn reflects the cosmological expansion rate at that epoch. 

In terms of the dark energy properties, the three points above mean that the dark energy equation of state is probed 
at all the relevant redshifts where it significantly affects the cosmological expansion rate, but independently on its 
value today; this is a very important feature since it tends to remove the common problem which arises when trying 
to constrain dark energy models with the CMB, i.e. that the relevant effect is made at the lowest redshifts where 
all the models converge to the same behavior. We saw an example of such occurrence in the factor D above, where 
the distances in dark energy cosmologies are essentially dominated by the present dark energy equation of state value 
wq, which greatly reduces the sensitivity to the w{z) behavior at higher redshifts. Indeed, that is the reason why the 
CMB power spectrum is quite sensible to the effective dark energy equation of state Wef / , but poorly on its redshift 
evolution [s^- In the bispectrum, on the other hand, the signal comes from redshifts strictly greater than 1, cutting 
out the present, and this is certainly important to discriminate between the different models proposed to explain the 
evidence for cosmic acceleration. Specifically, if w{z) is different as in the models exposed in Section ^ we expect 
to see a difference in the bispectrum signal, greatly enhanced with respect to the one in the ordinary CMB power 
spectrum, as we show explicitly in the next Section. 

Before concluding, we comment on the role of the description we adopt for the non-linear tail of the density 
fluctuation power spectrum. In Fig. 2|we plot the ratio R between the non-linear and linear power spectra according 
to the recipe adopted here j4Qj, for the dark energy models we consider here. The three models yield a similar behavior 
in R, with some difference in the scale at which R rises for the onset of non-linearity, which in SUGRA occurs at 
lower wavenumbers for a fixed z; as we see in the next Section, that is not induced by the particular recipe adopted 
for the non-linear mapping of the spectrum, but merely reflects the linear perturbation dynamics, inducing at any 
given epoch a larger fluctuation amplitude in SUGRA with respect to the RP and DEg// cases, slightly shifting the 
non-linear rise of the spectrum toward smaller multipoles. On the other hand, we will see that the mechanism leading 
to the sensitivity of the bispectrum to the dark energy equation of state redshift behavior is different, depending 
substantially only on the time derivative of the linear perturbation spectrum. 



In this Section we illustrate the phenomenology described above in the dark energy scenarios defined in Section ITll 
by computing explicitly the bispectrum signal as a function of the multipole I. On small multipoles, i.e. large angles, 
the corresponding scales are outside the horizon, in linear regime, and Q{1) > 0. On the other hand, on sufficiently 
large multipoles, Q{1) probes sub-horizon scales where the non-linear regime dominates, and it has negative sign. 
The transition is located at the angular scale where the contributions from the scales in linear and non-linear regime 
balance in the integrand, so that Q (/) is zero. For simplicity, here we evaluate the bispectrum by taking equal 
multipoles, li = I2 = I3 — I- This is equivalent to probe the CMB pattern with equilateral triangles of different sizes 
corresponding roughly to 180// degrees; that is certainly not the only choice, but it is enough for our purposes here. 
Eq. (fT7|) simplifies as 



I 



kNL {zbs) , 



(20) 



r{zbs) 



V. SOLVING FOR w{z) WITH THE CMB SPECTRUM AND BISPECTRUM 




(21) 



The bispectra in the three dark energy models described in Section ^ are plotted in Fig. El In Fig. El we plot the 
corresponding CMB angular power spectra, together with the binned WMAP data. As it is evident, the models are 
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almost degenerate in the CMB spectra, but can be resolved by looking at the bispectrum. Indeed, for the reason we 
explain in detail below, the angular location at which the linear and non-linear bispectrum contributions balance in 
the SUGRA model is shifted by several tens of multipoles toward smaller angular scales with respect to the RP and 
DEeff cases, while in the ordinary CMB spectrum the shift is of the order of a few multipoles; this indicates that the 
bispectrum is almost 10 times more effective than the ordinary CMB angular power spectrum in discriminating the 
redshift dependence of the dark energy equation of state. 

The three models yield equal expansion rate at the present, but as we pointed out in the previous Section, the 
bispectrum reflects the perturbation dynamics, and thus the cosmological expansion rate, in the interval 0.1 ;$ z < 2, 

independently on the present, in the multipole range 1000 <L I 100. In Fig. [7|we show the perturbation growth 

factor H13|) in the three models, normalized to the same value at present. As it is evident, the curve corresponding to 
the SUGRA case is sensibly larger at almost all epochs. Correspondingly, the change at low redshift, i.e. when the 
dark energy starts to dominate the expansion, is stronger in the SUGRA scenario, inducing a larger amplitude for 
the dP^s, /dz > tail in the line of sight integral H18(l . For the same reason, the bispectrum contribution coming from 
the gravitational bispectrum time derivative on non-linear scales has lower amplitude, simply because the rise due to 
the onset of non-linearity has to overcome the gravitational potential decay, which is stronger in the SUGRA with 
respect to the RP and DE^ff cases. 

The net effect is that, for a given multipole /, the larger is the value of w at the relevant epoch, the larger is the 
contribution from the linear regime decay of the gravitational potential. Thus, the scale at which the non-linear 
power balances the linear one is shifted toward larger wavenumber, i.e. larger multipoles, as in Fig. [S] Specifically, 
the cusp is located at I = 361,362 for the DEe// and RP models, respectively, and at I — 412 for SUGRA. As we 
already mentioned, the corresponding shift in the CMB angular power spectrum in FigEl is just a few multipoles. 
Note that even if we increase the present equation of state in the RP case to the border of the current constraints, 
say at wq = —0.78, the cusp moves to / = 390, still well below the SUGRA case, simply because no RP Quintessence 
model, within the current constraint, yield a w{z) behavior as the SUGRA one in Fig. ^ 

We conclude that the effect we pointed out in the previous Section represents a new and relevant feature induced by 
the dark energy dynamics on the CMB, which is worth taking in serious consideration in future observations: it causes 
a shift of several tens of multipoles in the CMB bispectrum arising in Quintessence scenarios yielding at the present 
the same equation of state, the latter being well within the range allowed by the current constraints. The strength of 
this result is evident when comparing the CMB spectra and bispectra in the different scenarios; the projection shift 
in the bispectra is about a factor 10 stronger than the corresponding one in the angular shift in the ordinary CMB 
spectrum. 



VI. SUMMARY AND CONCLUDING REMARKS 



The recent evidence for cosmic acceleration is one of the greatest surprises in modern cosmology, and imposes 
the introduction of a gravitationally repulsive force acting on cosmological scales, known as dark energy 1, 3]. In 
the wide set of ideas proposed to explain this occurrence (see [T^ ) , a class of models describes the dark energy as a 
self-interacting scalar field, the Quintessence, featuring a potential which provides the vacuum energy density required 
to explain the observed level of cosmic acceleration; several shapes for the Quintessence potential have been proposed 
out of super-symmetry/gravity breaking models |15lll6l|. These scenarios predict characteristic time variations of the 
dark energy equation of state w . IS j . By getting data on the redshift evolution of w throughout the relevant interval 
for cosmic acceleration, say z ^ 2, we can greatly restrict the range of models allowed to explain the dark energy and 

cosmic acceleration, getting constraints on the fundamental physics behind this process [l9l | . 

From this point of view, it is important to focus on the existing cosmological observables and build up new ones which 
are able to probe the cosmological expansion during the onset of cosmic acceleration. The obvious, cosmological process 
which occurs at similar redshift is the structure formation, that can be probed statistically by the weak gravitational 
lensing distortions which makes on background light 0, [s^. The latter effect has been recently observed and 
represents a promising probe for cosmology [23, |3| • In this work we focus on the weak lensing effect induced by 
structure formation on the background CMB anisotropics [25l |2^ . We look at the bispectrum, built out of the third 
order moment of CMB anisotropy statistics, because if the latter is Gaussian at last scattering as it appears from 
recent CMB observations p6| . then the third moment is zero within cosmic variance, with additional power imprinted 
by the weak lensing of the structures along the line of sight of CMB photons: the latter is expressed as a redshift 
integral between the last scattering and the present involving cosmological distances and the redshift derivative of the 
cosmological gravitational potential, induced the structures themselves |25l 123. 129L IsoL l4^ . 

We study the redshift behavior of the bispectrum signal arising from the correlation between weak lensing and 
ISW effect; that is made by two main contributions, deriving from the linear and non-linear regimes of the density 
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FIG. 5: Absolute value of dimensionless bispectrum for DEg// (solid line), RP (dotted line) and SUGRA (dashed line) with 
Wo = —0.9 as function of multipole /. 




FIG. 6: Dimensionless power spectra for our three models as function of multipole I, notice that the degeneration between 
models is greater than the bispectrum. All the models are normalized to WMAP data. Error bars are the binned WMAP data. 
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FIG. 7: Growth factor as function of redshift for our three models: DEe// (solid line), RP (dotted line) and SUGRA (dashed 
line). We have an appreciable difference between SUGRA and other two models in spite of the same value of wo. 



perturbations, yielding opposite signs after the horizon crossing: while the hnear regime makes the gravitational 
potential decreasing in time, the non-linear one, which can be described accordingly to the existing semi-analytical 
recipes mEll, yields a growth. Moreover, the integrand gets to zero at high and low redshift, and the relevant 
contribution, accounting for both the linear and non-linear regimes, comes from a narrow redshift interval centered 
around the typical epochs of structure formation, coinciding with the onset of cosmic acceleration. As a result, the 
CMB bispectrum is quite sensitive to the cosmic expansion rate at that redshift, i.e. on the dark energy equation of 
state at that time, independently on its present value. That represents a new and important aspect of the CMB as a 
dark energy probe, going in the direction of removing the effect coming from the dark energy evolution at the lowest 
redshifts, where all the models have to converge in order to reproduce the observed amount of cosmic acceleration 
|33j| . To illustrate this effect, we considered two Quintessence scenarios both having present equation of state close 
to —1, but markely different in the past, as a result of the different potential shapes: the shallow inverse power law 
and the steep one including an exponential correction. The latter scenario generally predicts an abrupt change of 
the dark energy equation of state at z ;S 1, dropping from values which can be as high as —0.2, while in the first one 

the w dynamics yield variations at the 10% level. As a result, the gravitational potential decay in linear regime gets 
enhanced in the SUGRA model with respect to the RP one; correspondingly, the scale at which the non-linear part 
of the spectrum balances the linear one is shifted toward gets smaller, subtending smaller angles in the sky. This 
induces a projection feature making the bispectra in the two scenarios shifted by several tens of multipoles, i.e. 10 
times stronger than the corresponding shift in the ordinary CMB power spectrum. 

We did not address here how much this effect is robust against variations in the other cosmological parameters. 
However, since it is mainly a projection feature, we argue that it could be mimicked by variations in the cosmological 
parameters inducing geometric features in the CMB anisotropics, such as the overall cosmic geometry, Hubble constant 
etc. On the other hand, it is important to note that such variations would affect similarly both the angular power 
spectrum and the bispectrum. We remark that the effect we report here affects the CMB bispectrum much more 
efficiently than the CMB spectrum, and this is made possible by three key features which are difficult to reproduce 
by varying the other cosmological parameters: (i) the fact that the bispectrum signal induced by forming structures 
through weak lensing probes a narrow redshift interval which docs not include the present, (m) the unique prediction 
of the dark energy cosmology, i.e. that one of the cosmological parameter, i.e. the equation of state, is time varying, 
and (in) that different dark energy scenarios generally predict different time variations of the equation of state. For 
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this reason we think that it will be useful to take into account the effect we pointed out here when dealing with data 
provided by the next, high resolution and high sensitivity cosmological probes such as Planck, CMBpol, as well as 
the weak lensing surveys which will be carried out in the next decade. 
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APPENDIX A: WIGNER 3J SYMBOLS 

The Wigner 3J symbols (W3J) come out from the composition of three spherical harmonics 



{2h + 1) {2I2 + 1) (2^3 + 1) f h h h \ ( h k h 



4tt \ / \ nil m2 



(Al) 



Eq. lAip represent the Gaunt integral. The W3J has a closed formula only for a few combinations of li, I2 and I3 and 
TOi. 7712 and ma; for the most general case we need a recursion relation, to be able to compute the whole set of W3J. 
The recursive relation for W3J is computed for fixed li, k, rrii and for the other possible values of 7712; m^is given 
by the relation 7713 = — (mi + 1712) and if this equality is not verified, the W3J are identically zero. The recursive 
relation is 

C (ma + 1) g (ma + 1) + D (ma) g (ma) + C (ma) 5 (ma - 1) = , (A2) 

where 

9M=( ^'), (A3) 

are the W3J as function of ma, and the coefficients C(ma) and Z?(ma) are defined by 



C(ma) = ^J{l2 - rna + 1) {h + ma) (^3 + m3 + 1) (^3 - m3), (A4) 
D(m2) = /a (/a + 1) + ^3 (^3 + 1) - ^1 ('i + 1) + 2m2m3. (A5) 

The usual choice for the starting values of W3J at the boundaries is g (mamin) = 9 (ma max) = 1 where mamin and 
ma max are respectively the smaller and the largest values of ma , are 

mamin = max[-;a,- (?3 +mi)] , (A6) 
ma max = min [^3,^3 - mi] . (A7) 

It is also possible to work only with the forward recursion, in this way it is sufficient to normalize the whole set of 
symbols with the normalization condition 



m2 mas 



_ (2?i + l).g^(m0 - 1 , (A8) 

rrii— 7712 min 

and the phase convention 

sign [g (ma max)] = (-1)'^-'^'-'"^ . (A9) 

3 

The W3J have an analytical expression in terms of factorials, when mi — ma = m3 — and L = ^ li is even: 



i=l 



h h Z3\ (V2)! A yiT^M 
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For odd L the WSJ is zero. Since the computation of factorial may be hard at high multipoles, we can use the Gosper 
factorial approximation: 



2n+\\^(-Y . (All) 



3J Ve 

in this way we can rewrite Eq. (|A10|I in a computational simpler form 

h h l3\ ( L \^''^ 1 /3e3L + l (6L-12;, + 1) 



0/ \ L + l) (6L + 7)i/n^ ^+1 / (3L-6Z, + 1)1/2 
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